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We investigate the hypothetical process of gravitational Cherenkov radiation, which may occur 
in modihed gravity theories. We obtain a useful constraint on a modihed dispersion relation for 
propagating modes of gravitational waves, which could be predicted as a consequence of violation of 
the Lorentz invariance in modified theories of gravity. The constraint from gravitational Cherenkov 
radiation and that from direct measurements of the gravitational waves emitted by a compact binary 
system are complementary to each other. 


Gravitational waves, predicted in gravitational theo¬ 
ries, are at the frontier of physics and astronomy. Re¬ 
cently, a variety of the gravitational theories have been 
proposed, motivated by a possible explanation of the ac¬ 
celerated expansion of the universe. For example, in 
recent research viable gravitational theories of gravita¬ 
tional waves with a mass term have been proposed [1-3]. 
Furthermore, graviton oscillations are predicted in the 
ghost-free bigravity theory [4, 5]. In the most general 
scalar tensor theory of the second derivative [6], which 
was rediscovered recently [7], the propagation speed of 
gravitational waves may deviate from the propagation 
speed of light [8]. Within general relativity, the propa¬ 
gation speed of gravitational waves is the same as that 
of light, whose deviation is related to the breaking of 
Lorentz invariance [9, 10]. Thus, the gravitational wave 
is important to characterize modified theories of gravity. 

The dispersion relation for gravitational waves has 
been much argued in the literature. For example, con¬ 
straints on the mass of gravitational waves have been dis¬ 
cussed [11]. Although gravitational waves have not been 
directly detected, the progress of observations such as 
those from the advanced LIGO project and the KAGRA 
project will make direct measurements possible. Assum¬ 
ing such future prospects of observational experiments, 
Mirshekari, Yunes, and Will investigated a future con¬ 
straint on the modified dispersion relation through grav¬ 
itational wave experiments [12]. They demonstrated that 
a stringent constraint on the mass term can be obtained 
with future direct measurements of gravitational waves 
emitted by a compact binary system. The authors of 
the recent papers [13, 14] further investigated the orbital 
evolution of binary pulsars in modified gravity models 
with Lorentz symmetry breaking and obtained a useful 
constraint on the model parameters. 

However, it is known that gravitational Gherenkov ra¬ 
diation (GGR) arises if a massive particle moves faster 
than the speed of the gravitational waves, which is pos¬ 
sible when the propagation speed of the gravitational 
waves is smaller than that of light [15]. This hypotheti¬ 
cal process puts a stringent constraint on the propagation 
speed of gravitational waves by including observations of 


extremely high energy cosmic rays [16]. The usefulness is 
indeed demonstrated for modified gravity models of the 
new Ether-Einstein gravity and the Galileon-type cosmo¬ 
logical model [17, 18]. 

In the present letter, we investigate the constraint on 
the modified dispersion relation for gravitational waves 
from GGR. We assume the same modified dispersion re¬ 
lation for the gravitational waves in Ref. [12]: 

^k = + Ak‘^U^, ( 1 ) 

where fc, Cg, and trig are the wave number, the prop¬ 
agation speed, and the mass of the gravitational wave, 
respectively, and A and a are also the parameters of the 
modihed dispersion relation. In the absence of the terms 
of TTig, A, and a, it is known that GGR puts the following 
constraint on the propagation speed [16]: 
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However, the GGR process for the modihed dispersion re¬ 
lation of the form (1) has not been discussed. We demon¬ 
strate that GGR puts a stringent constraint on A and Cg, 
depending on a and the sign of A, and that GGR cannot 
put a useful constraint on rrig. From the direct mea¬ 
surement of gravitational waves, the constraint on rrig is 
stringent, but the constraint on A is not very stringent 
[12], which means that the two methods are complemen¬ 
tary. 

Following Ref. [16, 18], the rate of emitting energy 
through GGR from a massive particle with mass m and 
relativistic momentum p = |p| can be written as 

= GnP^ / dkk^ / dcosffsin"^^ 
at Jo 7-1 

X Hy tOk'), (3) 

where Gn is Newton’s universal gravitational constant, 
Uli and ilf are dehned by + "m? and 

= \/(p-k)2 -\-nn?, respectively, and 5d is the 
Dirac delta function. Using the identity S]j{Uli — Ulf — 
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FIG. 1. Constraint on the model parameters from GCR for 
mgjp and 5 in the case A — 0. The shaded region satisfies 
(11); this region is excluded by GCR, where we assumed p = 
10^^ GeV and ct = 1 Mpc. The dashed curve shows (8) 
with k — p, above which the GCR process does not appear 
effectively. 


ujk) = 20/i5_D(r2j — {fli — uikY), we may write 


where we set = P- We estimate the condition at 
which the damping from GCR is significant for an ex¬ 
tremely high energy cosmic ray with initial energy p dur¬ 
ing time t as dE/dt > p/t, which yields 
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Integration of (3) makes a nontrivial contribution if 
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which is a condition at which GCR arises. Assuming 
m/p < 1, mg/k < 1, A(csfc)““^ < 1, and |(5| <C 1, 
where we defined <5 = 1 — c^, we have 


The shaded region in Fig. 1 satisfies this condition (11), 
which is excluded as the constraint of GCR. Here we 
assume an extremely high energy cosmic ray of p = 
10^^ GeV from a cosmological distance ct = \ Mpc. 
The dominant contribution of the k integration in Eq. (9) 
comes from the region of k of the order of p. The dashed 
curve in Fig. 1 shows the border of (8) with k replaced 
by p, i.e., 6 = m^/p'^. In the region above the dashed 
curve, GCR does not occur effectively. Thus we obtain 
the constraint from GCR process for S ^ 10“^®; however, 
the constraint on mg is of the order of mg ;$ 10^ GeV. 
When a high energy cosmic ray of p = 10^ GeV from a 
cosmological distance ct = IMpc is assumed, we obtain 
the constraint on the mg is of the order of mg Y 10 GeV 
for 6 Y 10“^. Thus the constraint on the graviton mass 
from GCR is not very stringent. 

Next we consider the case mg = 0; the condition at 
which GCR arises, (7), is 
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— p < 0, and then condition (5) leads 
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First we consider the case A = 0 but with finite gravi¬ 
ton mass mg Y 0 Cs < 1, i.e., (5 > 0. In this case, (7) 
gives 





( 8 ) 


(5-Afc“-^>0. (12) 


One can observe that GCR arises even when (5 < 0, i.e., 
Cs > 1, if A is negative. This is quite a contrast to the 
usual conclusion that GCR arises only when < 1 [15, 
16]. We integrate Eq. (3), adopting the approximation 
0 <C 1, which gives 
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thus GCR does not appear when ^ = 0. We integrate We here estimate the condition at which no significant 
Eq. (3) by adopting the approximation 0 <C 1, and we radiation damping occurs from GCR for an extremely 
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FIG. 2. Constraint on the model parameters from GCR for 
the case a = 3 and rrig = 0. The oblique straight line is 
the border of (12) with k replaced by p, and GCR does not 
occur in the upper region effectively. The elliptic curve is 
the border of (15), inside of which (dark-gray region) satisfies 
(15), where the GCR effect is so small that the extremely 
high energy cosmic rays do not damp. The light-gray region 
is excluded by GCR. 



FIG. 3. Same as Fig. 2, but for the case a = 4. 


high energy cosmic ray with initial energy p during time 
t as 


dE ^ p 
dt ~ t' 


(14) 


which gives 
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We examine conditions (12) and (15). Here we as¬ 
sume that the extremely high energy cosmic rays of 
p = 10^^ GeV from a cosmological distance ct = 1 Mpc 
are not signihcantly affected by the damping through 
GCR [16]. Figures 2 and 3 show the constraint on the 
model parameters from GCR for cases a = 3 and 4, re¬ 
spectively. The oblique straight dashed line in each figure 
shows the border of condition ( 12 ) with k replaced by p, 
above which GCR does not arise effectively. GCR arises 
in the lower region of this straight line. The elliptic re¬ 
gions in Figs. 2 and 3 satisfy condition (15), where GCR 


arises but the effect is small. The lower light-gray region 
in each figure is excluded by the GCR process. 

For the case (5 = 0, i.e., Cs = 1, (15) yields 
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Thus, for the case Cs = 1, by combining (12) and (15), 
the constraint from GCR excludes the value 
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for a = 4, respectively. 

Note that GCR excludes the negative region of the 
parameter of A when (5 = 0, i.e., Cs = 1. When (5 0, 

the constraint on A depends on 6. When (5 > 0, a positive 
region of A is excluded depending on the value of 6. 

Now we discuss the implication of our results for spe¬ 
cific models with the modified dispersion relation. First, 
we consider the model in broken-symmetry scenarios 
[19, 20], in which Lorentz invariance is broken associated 
with the Planck energy scale: 
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where Ep is the Planck energy scale and 77 is a nondimen- 
sional parameter. This model corresponds to our model 
parameters a = 3 and A = p/Ep and (5 = 0. From (17), 
the following region is excluded: 
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The authors of Ref. [19] argue that p naturally takes a 
negative value in the broken-symmetry scenarios, and so 
GCR puts a useful constraint on this scenario. In the 
previous work in Ref. [21], a constraint similar to ours 
is obtained for the photon’s modified dispersion relation, 
by comparing synchrotron radiation in the presence of 
the modified dispersion relation and observations from 
the Crab Nebula. 

We next consider the scenario with an extra dimension 
[22], which is an example of the case oi a = 4: 

a;2 = A:2+m2-aedt||, (21) 

where oiedt is a nondimensional parameter. This model 
corresponds to our model parameters a = 4 and A = 
—Oiedt/Ep and (5 = 0. The constraint from GCR (18) 
excludes the value 
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GCR constrains a positive value of aedt; otherwise, GCR 
does not occur. In this scenario of the modified dispersion 
relation, the sigir of aedt is not specified. The authors of 
Ref. [22] argue that the case aedt < 0 is relevant to theo¬ 
ries in which a generalized uncertainty principle is used. 
Therefore, unfortunately, our results cannot resolve this 
argument. 

As discussed in Ref. [12], there are other possible mod¬ 
els that predict the modified dispersion relation. For ex¬ 
ample, the Horava-Lifshitz gravity model predicts a mod¬ 
ified dispersion relation [23-27]: 
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for the left-right polarization modes, where k, w, and 
rj are parameters. The condition (7) at which GGR arises 
can be generalized as follows: 
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(24) 


Following the theory with the detailed-balance condition, 
ry = 0, the condition (24) is not satisfied as long as d < 
0, which does not allow GGR process. For the theory 
without the detailed-balance with 77 > 0, the condition 
(24) might be satisfied even when 6 = 0. For simplicity, 
we consider the case /r = 0, which reduces to an example 
of the case a = 6 with A = k^/4w^ — k?r\l2. From (15), 
we have the constraint on 6 and A. For the case 5 = 0, 
GCR excludes the value A{= hA — i<Ari/2) < —7 x 
10“®°GeV“^, by adopting p = 10^^ GeV and ct = I Mpc. 

The theory with noncommutative geometries, dis¬ 
cussed in Ref. [28], predicts a modified dispersion rela¬ 
tion of the form (1) with a = 4 at the lowest order of 
perturbative expansion of the momentum. The sign of 
the parameter A is naturally predicted to be positive. In 
this case, GCR does not appear as long as Cg = 1. 


In Ref. [12], future prospects of constraining the modi¬ 
fied dispersion relation from direct measurements of grav¬ 
itational waves from binary systems with the Advanced 
LIGO Project are addressed. With the method proposed, 
it is possible to put a stringent constraint on the graviton 
mass rUg. For example, for the case a = 3, the authors 
of Ref. [12] predict that the constraints at a level of 
mg < 10“^^-10“^^ eV will be obtained, depending on 
the target and the signal-to-noise ratio. However, their 
constraint on A is not stringent. The best constraint is 
JA] < 10 GeV“^ for the case a = 3. However, the cor¬ 
responding constraint on A from GCR is quite stringent, 
A A — 10~^'^ GeV~^, although the constraint on mg is 
not very stringent. Thus, the two methods are comple¬ 
mentary. 

We have investigated the constraint on the modified 
dispersion relation for gravitational waves from GCR, as¬ 
suming the form (1). The constraint on mg is not very 
tight; however, the constraint on A from GCR can be 
very stringent. When the propagation speed is less than 
the velocity of light, Cg < 1, GCR may appear even when 
A is positive, which is constrained through the GCR pro¬ 
cess. When Cg = 1, GCR puts a constraint on A only 
when A is negative. This constraint is very stringent 
compared with that from direct measurements of gravi¬ 
tational waves from binary systems as found in Ref. [12] 
(see also [13, 14]). Thus, the constraint from GCR is com¬ 
plementary to that from direct measurements of gravita¬ 
tional waves. 

After this paper was nearly completed, Ref. [29], in 
which a similar problem is investigated, appeared on the 
arXiv. We focused on the modified dispersion relation 
of the form (1) but the authors of Ref. [29] investigate 
constraints on various Lorentz violation operators. 
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